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ABSTRACT 

The angular power spectrum estimator developed by Peebles (1973) and Hauser 
\^ • & Peebles (1973) has been modified and applied to the 2 year maps produced by the 

COBE DMR. The power spectrum of the real sky has been compared to the power 
psj . spectra of a large number of simulated random skies produced with noise equal to the 

^ \ observed noise and primordial density fluctuation power spectra of power law form, 

with P{k) oc A;". Within the limited range of spatial scales covered by the COBET^MR, 
corresponding to spherical harmonic indices 3 < £^30, the best fitting value of the 
f^ ' spectral index is n = 1.2^~^q\^ with the Harrison-Zeldovich value n = 1 approximately 

"^ ■ 0.5o- below the best fit. For 3 < £<19, the best fit is n = 1.46+o:44- Comparing the 

^ '- COBE DMR AT/T at small (. to the AT/T at ^ « 50 from degree scale anisotropy 

o . experiments gives a smaller range of acceptable spectral indices which includes n = 1. 

o 



'i 



1. Introduction 

The spatial power spectrum of primordial density perturbations, P{k) where k is the spatial 
wavenumber, is a powerful tool in the analysis of the large scale structure in the Universe. In the 
first moments after the Big Bang, the horizon scale ct corresponds to a current scale that is much 
smaller than galaxies, so the assumption of a scale free form for P{k) is natural, which implies a 
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power law P{k) oc /c". Harrison (1970), Zeldovich (1972), and Peebles &: Yu (1970) all pointed out 
that the absence of tiny black holes implies n ^ 1 , while the large-scale homogeneity implied by 
the near isotropy of the Cosmic Microwave Background Radiation (CMBR) requires ntl. Thus 
the prediction of a Harrison-Zeldovich or n = 1 form for P{k) by an analysis that excludes all 
other possibilities is an old one. This particular scale-free power law is scale-invariant because 
the perturbations in the metric (or gravitational potential) are independent of the scale. The 
inflationary scenario of Guth (1981) proposes a tremendous expansion of the Universe (by a factor 
> lO^'^) during the inflationary epoch, which can convert quantum mechanical fluctuations on a 
microscopic scale during the inflationary epoch into Gpc-scale structure now. To the extent that 
conditions were relatively stable during the small part of the inflationary epoch which produced the 
Mpc to Gpc structures we now study, an almost scale-invariant spectrum is produced (Bardeen, 
Steinhardt Sz Turner 1983). Bond & Efstathiou (1987) show that the expected variance of the 
coefficients aim in ^ spherical harmonic expansion of the CMBR temperature given a power law 
power spectrum P{k) oc A;" is < a|„ > ocT[i+{n- l)/2]/r[^ + (5 - n)/2] for £ < 40. Thus a study 
of the angular power spectrum of the CMBR can be used to place limits on the spectral index n 
and test the inflationary prediction of a spectrum close to the Harrison-Zeldovich spectrum with 
n = 1. 

The angular power spectrum contains the same information as the angular correlation function, 
but in a form that simplifles the visualization of fits for the spectral index n. Furthermore, the 
off-diagonal elements of the covariance matrix have a smaller effect for the power spectrum than 
for the correlation function. However, with partial sky coverage the multipole estimates in the 
power spectrum are correlated, and this covariance must be considered when analyzing either the 
correlation function or the power spectrum. 

The power spectrum of a function mapped over the entire sphere can be derived easily from 
its expansion into spherical harmonics, but for a function known only over part of the sphere this 
procedure fails. Wright (1993) has modified a power spectral estimator from Peebles (1973) and 
Hauser &: Peebles (1973) that allows for partial coverage and applied this estimator to the DMR 
maps of CMBR anisotropy. We report here on the application of these statistics to the DMR 
maps based on the first two years of data (Bennett et al. 1994). Monte Carlo runs have been used 
to calculate the mean and covariance of the power spectrum. Fits to estimate {Q\ms)^^ ^^^ "■ 
by maximizing the Gaussian approximation to the likelihood of the angular power spectrum are 
discussed in this paper. Since we only consider power law power spectrum fits in this paper, we 
use Q as a shorthand for {QJims)^'^ °^ Qrms-ps, which is the RMS quadrupole averaged over the 
whole Universe, based on a power law fit to many multipoles. {Q\ms)^^ should not be confused 
with the actual quadrupole of the high galactic latitude part of the sky observed from the Sun's 
location within the Universe, which is the Qrms discussed by Bennett et al. (1992a). 



2. Estimating the Angular Power Spectrum 



Wright (1993) has discussed the modification of the Hauser-Peebles angular power spectrum 
estimator for use on CMBR anisotropy maps. We include a description of this method for 
completeness. Consider a collection of spectral functions F^^ which are defined to be orthonormal 
in the measure dQ/Air. These are the real spherical harmonics, normalized to have an RMS value 
of unity for each harmonic. 

The inner product of spatial functions / and g is defined as 

<fg>= J f{n)g{n)dn/47r. (i) 

Note that the F^m's satisfy 

< FemFe'm' >= ^e ^m ■ (2) 

Given the temperature distribution T(r2), define the RMS power at each multipole i as 

e 
Tf= J2 \<FemT>\^ (3) 

m=—£ 

The Hauser-Peebles approach to power spectra on the sphere with non-uniform or absent coverage 
involved correcting for the average density of sources: the Fqo term. In the case of the DMR maps, 
we clearly should also correct for dipole terms. 

Redefine the inner product to apply to the non-uniformly covered sphere: 



< h >= 5|P^ (4) 






where j is an index over pixels, and Wj is the weight per pixel. In the galactic plane, Wj = 0. The 
galactic plane cut used in this paper excludes the 1/3 of the sky with \b\ < 19.5°. Outside of the 
galactic plane, one can choose whether to have Wj follow the map weights based on the number 
of observations, N^bs- We have used uniform weights instead of N^bs weights, which increases the 
effect of radiometer noise in the results, but also reduces and simplifies the correlation between 
different T^'s in the result. 

Now define revised functions Gim for ^ > 1 given by 



-^ 00 < FpoFem > _ ■^ Fl,m' < Fl,m'Flm > 



Wm — i'im J-f=, — ^ ^ 2^ ^-f^ ^ ^ • [o) 



These functions are orthogonal to the monopole and dipole terms in the region covered by the map 
with the specified weights. But they are not orthogonal to each other, nor are they normalized. 
The function G31 is substantially affected by the dipole removal, since the galactic plane cut 
couples harmonics with Ai = ±2 and Am = 0. On the other hand Gi3^i3 is not much affected by 
the monopole plus dipole removal but is far from normalized in the polar caps, since most of its 
power is concentrated in the galactic plane. 



We can now define the terms used by Hauser & Peebles: the normahzation integral for G^m 

JP =< GirnGim > (6) 

and the estimated spectrum 

zr = \< GemT > iv^r- (7) 

Hauser & Peebles recommend that the estimate to be used for the spectrum should be the average 
value of the Z^^s weighted by the J^^'s. This quantity is 

'^' ./ 7m \ Zm=-<; < GjmT > 

'< ^i >m- =r« — • i»j 



The DMR experiment has two independent channels, Ta and Tg, at each of three frequencies: 
31, 53 and 90 GHz. The sum and difference maps formed from the A and B channel maps can 
be used to determine error associated with this estimate of T^^. Since T^^ is obtained as a sum of 
squares, it is necessarily positive, and is thus a biased estimator. The sum and difference maps 
can also be used to correct this bias. Let the sum map be 5 = (Ta + Tb)/2 while the difference 
\s D = {Ta — Tb)/2, where Ta and Tb are the maps produced by the A and B sides of the DMR 
instrument. Then an unbiased estimate of the true power spectrum of the sky is given by 

T-2 _ (^n I -,^llm=-l (< GimS > - < GgrnD > ) 

ii^[zi + i) —J — . (yj 

These statistics evaluated for the 53+90 GHz maps with the first 2 years of data are shown 
in Figure |l[ Assuming that both the noise map in D and the cosmic plus noise map in S are 
described by isotropic Gaussian random processes (independent of m), we get an estimate for the 
uncertainty in T^: 

, , 2i2i + ifj:{jrf\iT.<Gi^s>^f + {j:<GimD>Y 

^ i^l) = ^-—74 -■ (10) 



This error estimate provides the error bars in Figure |l|. The Monte Carlo simulations discussed 
below have shown that this error estimate is correct: the mean over many simulations of the 



variance in Equation 10 agrees with the variance computed from the scatter in the power spectra 
computed using Equation |9[ This uncertainty can easily be approximated for the case of no 
galactic plane cut, and small signal to noise ratio. In this case J™ = 1, and the expected value of 
< GimS >^ and < GimD >^ are both oc al/Ntot, where ui is the uncertainty in a single DMR 
observation and Ntot is the total number of observations over the whole sky. Thus the variance of 
T/ is 

a\T!)^{2i + l)^ (11) 

in this case. For the Harrison-Zeldovich spectrum predicted by inflation, the signal to noise ratio 
of T^ varies like i~^'^. Because of this rapid decrease of significance with increasing £, we have 



constructed binned statistics by summing neighboring T"^ into bins covering the ranges £ = 2, 
3, 4, 5-6, 7-9, 10-13, and 14-19. These bins are approximately uniform in ln£. These binned 
statistics are used on pfots to avoid clutter, but the maximum likelihood fits discussed below use 
the unbinned statistics. 

These Tf^s are quadratic statistics derived from the DMR maps. Wright et al. (1994) define 
an averaged response of a quadratic statistic to the spherical harmonics of a given order i' . Let 
'^fe.'m' b^ ^^^ response in the t order when the input is the spherical harmonic Fii^i. The mean 
over m' of this quantity, needed to analyze isotropic random processes, is 

Ti,i' - ^f^^ (12) 

Table |l] shows 1000 times this quantity for £ = 2 ... 19 and £' = . . . 19 when the galactic plane is 
cut at |6| = 19.5°. Note the strong coupling of orders separated by Ai = ±2 caused by the galactic 
plane cut. With no galactic plane cut, T^^, = 5^ for £ > 2. 

The values in Table [^ can be used to estimate the response to power law power spectra of 
primordial density perturbations with an amplitude Q and a power law index of n: 

J^TTrTZr. ., n2 y (2f + l) ^2 ^2 r[f + (n - l)/2]r[(9 - n)/2] 

T, {Q,n) ^ Q Z. ^^^^'^«' r[f + (5 - n)/2]r[(3 + n)/2] ^^^^ 

where Ge is the coefficient of the Legendre polynomial expansion of the beam given in Wright et 
al. (1994). The effective spherical harmonic index defined by Wright et al. (1994), 



^eff = 2 exp 



91n(r//Q2) 



dn 



(14) 



n=lJ 

can be evaluated either from the sum above or from the mean of Monte Carlo simulations. The 
result is that £(.ff is significantly smaller than £. The solid curve in Figure ^ shows the relationship 
for |6| > 19.5°. Even for the case of no galactic plane cut, i^jf is smaller than i when ^ > 2, as 
is shown by the dashed curve in Figure ||. For £'s beyond the DMR beam cutoff at £ ~ 19 the 
response to an n = 1 input spectrum is dominated by the off-diagonal response to low ts, so ^e// 
saturates. These high i statistics are primarily sensitive to high n models. 



Monte Carlo Simulations 



Monte Carlo simulations of the T^ statistics have been done for n = —0.75 to 2.75 in 
P{k) oc k"", and various values of Q. Since the power spectrum is a quadratic function of the sky 
temperatures, calculation at 3 different values of Q for a given realization of the detector noise 



and cosmic variance suffice to produce tlie result for all values of Q using quadratic interpolation. 
Therefore the power spectrum of a particular Monte Carlo realization is given by 

Tf{Q,n) = a{n)Q'^ + h{n)Q + c{n) (15) 



and the mean power spectrum of a set of Monte Carlo skies, T^{Q,n), is given by 



T^'{Q,n)=a{n)Q' + b{n)Q + c{n). (16) 



Note that the expected values of b{n) and c{n) are zero, but the actual values from a finite set of 
Monte Carlo simulations will be non-zero. The covariance matrix C{Q, n) of the Tf statistics is 
also determined using the Monte Carlo simulations. Since the Tf are quadratic functions of Q, 
the covariance matrix is a quartic polynomial in Q. The coefficients of the odd powers of Q in 
this polynomial have expected values of zero, so the covariance matrix breaks into a noise-noise 
part (the coefficient of Q^), a signal-noise part (the coefficient of Q^) and a signal-signal part (the 
coefficient of Q^)- Seljak & Bertschinger (1993) decompose the covariance matrix of the angular 
correlation function in the same way. 

The radiometer noise contribution to the simulated maps includes the positive noise 
correlation for pixels separated by 60° using a corrected version of the technique given in Wright 
et al. (1994). The DMR maps are found by solving the matrix equation AT = M (Lineweaver et 
al. 1994), where A is a sparse symmetric matrix, with diagonal elements An = Ni, the number 
of observations of the i'th pixel; and off-diagonal elements Aij equal to minus the number of 
times the i'th and j'th pixels were compared. Wright et al. (1994) assumed that the right-hand 
side vector M would be uncorrelated, but it is actually anti-correlated for pixels separated by 
60°. A correct way to generate correlated noise maps is to note that afA'^, with cji being the 
error in a single sample, is the covariance matrix of the noise maps. This implies that noise 
maps can be created using T = aiA^^'^U, where U is a vector of uncorrelated, zero mean unit 
variance Gaussian random numbers. Even though A is singular, a series expansion of A^^'^ 
converges rapidly except for the eigenvector corresponding to the mean of the map. This series is 
derived by writing A = D{I + E)D, with Dij = 6ij \/ An and Eij = (1 — bij)Aijl ^ AnAjj. Then 
^-0.5 r^ D^^(l — 0.5-E + CSTS-E^ — ...). The first term gives an uncorrelated noise map, while 
the second term gives a first-order correction for the 60° correlation that is exactly one-half the 
correction used by Wright et al. (1994). Thus we first generate a O'th order map using uncorrelated 
random numbers scaled by N^ ' . The first order correction is 1/2 of the weighted mean over the 
reference ring at 60° separation of the O'th order map values, with the weights given by the number 
of times each pixel pair is observed. The second order correction is 3/4 of the weighted mean 
over the reference ring of the first order correction. The m'th order correction is (2r7T, — l)/(2m) 
of the weighted mean over the reference ring of the [m — l)'th order correction. A similar series 
approximation for the covariance matrix itself is A^^ ~ D^^{I — E + E"^ — . . .)D^^. 



4. DMR Data Selection and Power Spectrum Estimates 

The data analyzed in this paper are the maps from the first 2 years of DMR data 
discussed by Bennett et al. (1994). The maps are made using pixels with cube faces oriented 
in galactic coordinates. To minimize the noise, a linear combination of the 53 GHz and 90 
GHz channels is made: O.6T53/O.93I + 0.4r9o/0.815. The denominators in this expression 
convert the Rayleigh-Jeans differential temperatures T53 and Tgo into thermodynamic AT's, and 
the 60:40 weighting is used because the 53 GHz channels are the most sensitive. This linear 
combination applied to the publicly released 1 year maps in ecliptic oriented pixels has also been 
analyzed. A cross-over version of this combination, using 53A+90B and 53B+90A, has also been 
analyzed. A second linear combination used is the "No Galaxy" map constructed using weights 
Ti\fG = — 0.4512r3i + 1.2737r53 + O.3125T90. This combination is calibrated in thermodynamic 
AT units, gives zero response to the mean galactic plane and to free-free emission (Bennett et 
al. 1992a). The 53 GHz maps are also analyzed by themselves, using T53/O.93I to convert to a 
thermodynamic AT scale. Finally, the cross-power spectrum of the 53 x 90 GHz maps has been 
found, by letting the sum map be the average of the 53 and 90 GHz maps, each converted into 
thermodynamic AT's, while the difference map D is (53 — 90)/2. Table ^ gives the binned power 
law statistics for these four data sets. The error bars are the square root of the diagonal elements 
of the binned covariance matrix C{Q,n) from the Monte Carlo runs, evaluated at the best fit 
values of Q and n, and thus include both radiometer noise and "cosmic variance" . The radiometer 
noise for each case is derived from the variance of the difference maps. The "cosmic variance" 
is the error in estimating the global mean properties of the Universe from a limited sample. It 



can be estimated from Equation IC and Equation |9| in the case where the difference map is zero, 
giving a limiting fractional precision of a{Tl)/Tf ^ y^ Air /['[}{£ + 0.5)], where J7 is the sky coverage 
(Scott, Srednicki & White 1994). 



Maximum Likelihood Estimation 



Given the mean power spectrum T^{Q, n), the covariance matrix C{Q, n) and the actual power 
spectrum Tf , define the deviation vector e^ = Tf — T^{Q, n) and the x^ statistic x^ = e^C~^e. All 
of the fits in this paper are based on the range i. = imin ■ ■ ■ ^max with £rnin = 3 and imax = 19or30. 
C is thus a 17 X 17 or 28 x 28 matrix. Ignoring the quadrupole is reasonable because the galactic 
corrections are largest for 1 = 2, and the maximum order used is set by the DMR beam-size of 
7° and the increased computer time required to analyze more orders. Since the magnitude of the 
covariance matrix gets larger rapidly when Q increases there is a bias toward large values of Q 
when minimizing x^- One can allow for this by minimizing — 21n(T) instead of x^, where L is the 
Gaussian approximation to the likelihood: 

- 2 ln(L) =x^ + ln(det(C)) + const. (17) 



Seljak & Bertschinger (1993) have applied this method to the correlation function of the DMR 
maps. This method has the interesting property that if the observed power spectrum matches the 
model exactly then the fitted value of Q is significantly less than the true value. At the minimum 
of x^) which is x^ = in this case, there is still a large slope in — 21n(L) because of the rapid 
variation of ln(det(C)) with Q. Figure ^ shows this effect: the diamond symbol shows the values 
of the amplitude and n obtained by minimizing — 21n(L) when the observed power spectrum is the 
mean of 4000 Monte Carlo's with n = 1 and Q = 17 fiK. It is clearly biased toward low amplitude 
when compared to the dots, which show fits to the individual power spectra from the 4000 Monte 
Carlo's. The big circle shows the results of minimizing — 21n(L) for the real sky power spectrum 
from the 53+90 map with 2 years of data. 

The maximum likelihood technique gives an asymptotically unbiased determination of the 
amplitude Q and index n, but only as the observed solid angle goes to infinity. Since we are 
limited to about 8 sr of sky, asymptotically unbiased means biased in practice. In addition, the use 
of a Gaussian approximation for the likelihood of our quadratic statistics can introduce additional 
errors. We use our Monte Carlo simulations to calibrate our statistical methods to avoid biased 
final answers. If we maximize L using the simplex method we find that the maximum likelihood 
index is biased upward from the input n used in the Monte Carlo's by ~ 0.1. An alternative 
method based on finding the zero in the finite difference L{Q, n + 0.2) — L{Q, n — 0.2) gave a much 
smaller bias but sometimes failed to converge for power spectra that were not well fit by a power 
law. 

The cross power spectrum of the real sky based on the (53A+90A) x (53B+90B) maps is best 
fit in the range 3 < ^ < 19 by an n = 1.55 model when we use the simplex method to maximize 
the likelihood. The improvement in 21n(L) between the fit with n forced to be 1 and the n = 1.55 
model is 2.5, which corresponds to 1.6 a. However, 14% of the Monte Carlo simulations made 
with riin = 1 and the maximum likelihood amplitude for n = \ give fitted values of n that are 
larger than 1.55, so this deviation from a Harrison-Zeldovich spectrum is really only a "1.06 
fj" deviation. Similarly, 54% of simulations made with nj„ = 1.5 and the maximum likelihood 
amplitude for n = 1.5 had fitted indices higher than 1.55, indicating that nj„ = 1.5 is actually too 
high by 0.10 o". Using the same procedure we find that nj„ = 0.5 is 2.07 cr low, nin = 2 is 1.41 o" 
high, and n^ = 2.5 is 2.81 a high. Interpolating to find values of n^ that deviate by -1, and +1 
a defines our quoted limits on the spectral index for 3 < £ < 19: n = lAQ'^^^^. 

With 4 years of data these limits will improve to An ~ _o 35 ^°^ the 53+90 maps if we assume 
that the maximum likelihood n remains the same. 

While waiting for this paper to be refereed, new computing facilities allowed us to increase 
^max to 30. The increased power at 20 < £ < 30 expected for Uin ~ 1.5 is not seen in the real 
maps, so the fitted values of n go down. Over the 3 < £ < 30 range the fits to the cross power 
spectra are n = 1.32+° j^ for (53A+90A)x(53B+90B), n = 1.22+°;^^ ^j. (53A+53B)x(90A+90B), 
and n = 1.20lo;46 for (53A+90B)x(53B+90A). These values have all been de-biased using the 



Monte Carlo simulations as discussed above. Each of these fits involves 4 of the 6 possible cross 
spectra among the 53A, 53B, 90A and 90B maps. Averaging these three cross spectra gives us all 
of the 6 possible cross spectra. The signal-to-noise ratio improvement from using 6 instead of 4 
cross products is quite modest, however, and is equivalent to a 22% increase in integration time. 
Thus the adopted range for the spectral index is n = 1.25lg;45. 

The spectral index from the NG maps is n = lAl^-^\j, where the large uncertainty is caused 
by the increased noise in the NG maps. The galaxy removal process subtracts the relatively noisy 
31 GHz channels of the DMR from a weighted sum of the quieter 53 and 90 GHz channels, and 
then rescales the result to allow for the partial cancellation of the cosmic AT by the subtraction. 
Both the subtraction and the rescaling increase the noise, and the overall process effectively 
doubles the radiometer noise. 

The difference between the n = 1.22_Q'4g reported here for 53 x 90 and the n = 1.15 reported 
by Smoot et al. (1992) is partly caused by the use of the real beam in this paper instead of the 
Gaussian beam approximation used by Smoot et al.. The ratio of G10/G4 from Wright et al. 
(1994) for the real beam to the same quantity for the Gaussian approximation is 0.92, and to 
compensate for the greater suppression of £ = 10 by the real beam the fitting procedure increases 
n by 0.2. This increase has been partly compensated by a decrease of n when going from the 1 
year to the 2 year maps. The de-biased fit to the 1 year (53A+90A)x(53B+90B) cross-power 
spectrum for 3 < ^ < 30 is n = imto.tt 

Bennett et al. (1994), using the real DMR beam instead of the Gaussian approximation, find 
that the maximum of the likelihood L{Q,n) occurs at Q = 12.4 ^K,n = 1.59 from an analysis 
of the cross-correlation function of the 2 year 53 x 90 GHz maps. This analysis included the 
quadrupole, and the low observed quadrupole leads to increased values of n when it is included in 
the fit. A no quadrupole fit gives the maximum likelihood at n = 1.21^055. 

Smoot et al. (1994) give estimates of the spectral index n derived from the variation with 
smoothing angle of the moments of the DMR maps, and of the genus of the DMR maps. The 
determination from moments is primarily based on the second moment, and the variation of the 
second moment with smoothing angle is equivalent to the power spectrum. This moment method 
gives n = l.Tlg'g when applied to the first year maps, which is quite consistent with the power 
spectrum of the first year maps. The genus method also gives n = 1.7 but n = 1 does not give a 
significantly worse fit. 

Gorski et al. (1994) examine linear statistics that are similar to (GimT). These have the 
major advantage that their distribution is exactly Gaussian, and thus the Gaussian form for the 
likelihood in Equation |l^ is exact. The linear statistics used by Gorski et al. define a position 
in a 961 dimensional space (for i < 30) which is hard to visualize, but using the exact Gaussian 
likelihood function for 3 < £ < 30, Gorski et al. (1994) find the maximum of L{Q,n) occurs at 
n = 1.02 for the combined 2 year 53 GHz plus 90 GHz map. Note that the 3 < ^ < 30 fits in 
this paper still include the small effect of the quadrupole on higher ts due to the off-diagonal 
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elements in the response matrix, while those in Gorski et al. (1994) are completely independent 
of the quadrupole. If Equation ^ is modified to also subtract quadrupole terms from the G^^'s, 
a different modified Hauser-Peebles power spectrum is obtained which is much more similar to 
the £ = 3 — 30 analysis of Gorski et al. (1994). In this variant the mean power in T4 for n = 1 
Monte Carlo skies goes down by 31% while T^ for the real sky goes up by 16%, leading to a higher 
1 = 4: point that balances the high £ = 14 — 19 bin and reproduces the Gorski et al. spectral index 
n = 1.0. Also note that Gorski et al. (1994) use a known shape and amplitude for the noise power, 
computed from the covariance matrix of the map, which allows them to use the "auto-power 
spectra" , while this paper just assumes that the A and B noises are uncorrelated and can only use 
cross-power spectra. 

Table ^ summarizes these results and includes model-dependent comparisons of COBE data 
to smaller-scale data. The apparent large-scale index we have used above is denoted Uapp, while 
the model-corrected primordial spectral index is n^^. 

Maximum likelihood fits for Q with n forced to be 1 over the 3 < ^ < 30 range give new 
determinations of the power spectrum amplitude: Q = 20.2 it 1.8, 19.2 it 1.3 and 15.6 it 2.1 fiK for 
the 53, 53+90 and NG maps. Comparing fits forced to n = 0.5 and n = 1.5 allows a determination 
of the effective wavenumber for these amplitudes: ieff = 6.1, 6.8 and 4.2 respectively. The 
amplitudes for the 53 and 53+90 maps are higher than the 17 /iK reported earlier because the 
maximum likelihood fit has emphasized the higher ts in determining the best fit since they have 
smaller cosmic variance, and shifting to higher ts gives a higher amplitude because the best fit 
value of n is greater than 1. The values from the NG maps are statistically consistent with the 
53+90 maps, but the possibility of a galactic contribution to Q and n is much reduced with the 
NG map. The maps are actually more similar than the 26% spread in best fit n = 1 amplitudes 
would suggest: the simpler o"(10°) statistic computed using DMRSMUTH (see Wright et al. 1994) 
in |6| > 30° is 31.9, 32.6 and 31.4 /xK for the 53, 53+90 and NG maps respectively, a spread of 
only 4%; while the GET_SKY_RMS program in |6| > 20° gives o-(10°) = 31.3, 29.1 and 30.7 /xK, a 
spread of only 7%. Thus most of the difference in the best fit amplitudes is caused by the shift of 
the weights to higher ts. 



6. Comparison ^vith Degree Scale Experiments 

Several groups have reported statistically significant signals from AT experiments with beam 
sizes and chopper throws close to 1°. These results are usually reported as limits on the amplitude 
of a Gaussian correlation function, 

Cg{e) = Cg{0) exp(-O.50V^c')- (18) 

We have calculated the conversion from the reported limits on Gaussian Cg{0) to limits on power 
law power spectra as follows: first, given the size of a Gaussian approximation to the experiment 
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beam, as =FWHM/V'8 In 2, find the beam smoothed Gaussian correlation function: 

Cg,sm{(^) = ggW ga+L^ e^P[-0-5^V(^c + 24)]- (19) 

Then the single subtracted, double subtracted, triple subtracted (Python) or square pattern (JI+) 
double subtracted (WD2) temperature difference is found from 

var(Ars5) = 2{Cg,sm{Q)-Cg,sm{e)) 

vaiiATDs) = l.5Cg,smiO) - 2Cg,,„(0) + 0.5Cg,s,ni2e) 

vaviATTs) = 1.25Cg,,„(0) - 1.875Cg,,„(0) + 0.75Cg,sm{29) - 0.125Cg,sm{3e) 

var{ATsQ) = Cg^smiO) - 2Cg,smiO) + Cg,smiV29) (20) 

where 6 is the chopper throw. The same temperature differences are then estimated for power law 
power spectra with Q = 17 /uK and n = 0.5, 1, and 1.5 using the expression for the beam smoothed 
power law correlation function 



^n,sm\y) — Q 



2 r[(9-n)/2] 
r[(3 + n)/2] 



With Q held fixed, we can define the effective spherical harmonic order for a given experiment 

using 

var(Ar[n = 1.51) , , 

where var(AT) is one of the four expressions in Equation |2^ with Cg^sm replaced by Cn,sm 
and the choice of which expression for var(AT) to use depends on the chopping strategy used 



in each experiment. The ratio of the observed var(AT) derived from Cg^sm in Equation 20 to 
var(Ar[n = 1]) with Q = 17 fiK. then defines the y-axis coordinate in Figure ^, while ieff gives 
the X-axis coordinate. 

Ganga et al. (1993, 1994) have have estimated the power spectral parameters Q and n using 
data from the Far Infra-Red Survey (FIRS), a balloon-borne survey experiment with a 3.8° beam. 
They obtain n = l.OTi'o which is consistent with the value in this paper. From the slope of 
their likelihood contours in the Q — n plane, we derive an £e// = 6.6 for their n = 1 amplitude 
determination, giving the FIRS point on Figure 0. 



7. Conclusions 

The Hauser-Peebles method of analyzing angular power spectra has been applied to the DMR 
maps. While the best fit to the observed power spectrum has n = 1.25Iq'45, ^^^^ deviation from 
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the n = 1 case is not statistically significant. To obtain a more accurate determination of n we 
need to compare the COBE DMR amplitude with ground-based and balloon-based experiments 
at smaller angular scales, which are sensitive to higher €s than COBE. 

The ULISSE and Tenerife experiments (Watson et al. 1992) with beam sizes near 6° and 
chopper throws of 6 — 8° give upper limits in the ieff ~ 15 range which support n < 1.5. The new 
Tenerife results (Hancock et al. 1994) at the same angular scale as Watson et al. give a central 
value that is slightly above the earlier upper limit. Both are plotted in Figure Q. 

The calculations of Kamionkowski &; Spergel (1994) suggest that for open Universes with 
il ~ 0.1 the power at low ^'s will be depressed relative to the n = 1 flat Universe prediction. This 
prediction is consistent with the data presented here, but fluctuations due to cosmic variance at 
low £'s are as large as the difference between the open Universe model and the scale-invariant flat 
Universe n = 1 model. 

The string model prediction given by Bennett, Stebbins & Boucher (1992b) also has lower 
power at small ^'s, and is thus consistent with the COBE angular power spectrum, but a cutoff at 
higher i is needed. Reionization of the Universe at redshift z will hide structures on scales smaller 
than 60°/\/l + z and provide the needed cutoff, but zt 100 is required to have a substantial optical 
depth with the baryon abundance derived from Big Bang nucleosynthesis (Walker et al. 1991). 
But reionization will not "smear out" the edges produced by strings seen at smaller redshifts. 
Thus, unlike most scientific models which can only be falsified, the string model can be verified by 
finding "the edge", which will remain infinitely sharp even with reionization. The sharp edges in 
the AT maps produced by nearby strings limits the slope of the cutoff to £~^ relative to an n = 1 
spectrum. Graham et al. (1994) find that the SP91 data is significantly non-Gaussian, which 
suggests that an edge may have been found. If true, this would increase the discrepancy among 
the degree-scale experiments, since the presence of an edge would increase the variance, but SP91 
has the smallest variance of the four degree-scale experiments. 

This expected increase in the variance due to non-Gaussian features is clearly present in the 
20 GHz OVRO and RING experiments which have the same angular scale. The RING experiment 
covered a larger region, and was contaminated by discrete sources whose existence was verified by 
the VLA. The 170 GHz MSAM experiment also saw what appeared to be discrete sources, and 
these were not included in the analysis. There is no sensitive, higher angular resolution telescope 
to verify that the large defiections seen by MSAM are indeed point sources. Thus it is possible 
that the large deflections are true cosmic AT's. The open circles above the MSAM data points in 
Figure |^ show the increased power that results if these data are not excluded in the analysis. 

The bulk flow data of Bertschinger et al. (1990) (at ieff ~ 10^) require n ~ 1 to agree with 
COBE (Wright et al. 1992), while the larger bulk flow on larger scales seen by Lauer & Postman 
(1992) requires n ~ 2.9 to agree with COBE, if we assume that the reported bulk flow represents 
the RMS velocity on this scale. 

The experiments at ~ 1° scale offer the possibility of a better determination of the primordial 
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power spectrum index n, but the model-dependent effects of the wing of the Doppler peak at 
i « 200 must be allowed for. Even in the large angle region ^ < 30 small model-dependent 
corrections must be made. In Figure |^, the upper Cold Dark Matter (CDM) curve has a primordial 
spectral index Upri = 0.96, but an apparent index Uapp = 1.1. Since the spectral index 1.25_q'45 
found in this paper is an apparent index, the COBE power spectrum is consistent with the 
prediction from inflation and CDM or Mixed Dark Matter models that Uapp ~ 1.1. On the other 
hand, vacuum-dominated models such as the Holtzman (1989) model with $7^ = 0.02, /i = 1, 
^CDM = 0.18, and ^vac = 0.8 will give Uapp = 0.9 for Up^i = 1.0 (Kofman & Starobinski 1985), 
which deviates by slightly less than \a from the COBE value. We have found the values of the 
primordial spectral index Upri that will connect the COBE NG amplitude and fg// found earlier 
with the degree-scale experiments using the scalar transfer function from Crittenden et al. (1993). 
We have ignored the tensor transfer function because the current accuracy in determining n is not 
sufficient to fix the tensor to scalar ratio, and because the excess quadrupole predicted by the 
tensor transfer function is not seen in the COBE power spectrum. The South Pole experiment of 
Schuster et al. (1993) at £ ~ 44 requires Upri « 0.48 it 0.34 to agree with COBE. The Saskatoon 
experiment of Wollack et al. (1993) at ^ ~ 55 requires Upn ~ 1.04 it 0.29 to agree with COBE. The 
PYTHON experiment of Dragovan et al. (1994) at £ ss 71 requires Upn « 1.58 ± 0.22 to match 
COBE. The ARGO experiment of de Bernardis et al. (1994) at £ « 75 requires np^i « 1.10 it 0.16 
to match COBE. The weighted mean of these values is Upri = 1.15 it 0.11. Unfortunately x^ = 8-0 
with 3 degrees of freedom when comparing these four values of Upri with this weighted mean, 
indicating that these four experiments are mutually inconsistent. If we allow for this discrepancy 
by scaling the error on Upri, we get a value Up^i = 1.15ib0.2ib0.1 from this comparison of COBE 
with the degree-scale experiments, where the second error bar is contribution of the uncertainty of 
the COBE NG amplitude to Upri. Thus this comparison of COBE with degree-scale experiments 
gives a more precise value the primordial spectral index that is still consistent with inflation. 
With more data from COBE (4 years are recorded) the large angular scale amplitude will become 
more and more certain. The £e// associated with this amplitude will shift to larger values ~ 10. 
Reliable, consistent determinations of 5T on scales £e// ^ 50 will be needed to compare with the 
large-scale AT. With only two years of data, the COBE DMR large scale amplitude has relative 
errors that are two times smaller than the errors of the current degree-scale experiments. Thus 
the degree-scale experiments need to be extended to a sky coverage that is 10 times higher than 
their current coverage to match the expected COBE uncertainty with four years of data, or else 
achieve an equivalent increase in accuracy by reduced noise or systematic errors. 

We gratefully acknowledge the many people who made this paper possible: the NASA Office 
of Space Sciences, the COBE flight operations team, and all of those who helped process and 
analyze the data. In particular we thank Tony Banday, Krys Gorski, Gary Hinshaw, Charlie 
Lineweaver, Mike Hauser, Mike Janssen, Steve Meyer and Rai Weiss for useful comments on the 
manuscript. 
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Table 1; Responses T^^(out) as a function of ^j„ normalized to 1000 
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Table 2: Ratio of the binned power spectrum from Equation 9 to a Q = 17 ^K, n = 1 model. 
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-351 
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229 
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Table 3: Binned Hauser-Peebles power spectrum of the 2 year NG maps, the best fit n = 1.4 
model, the nominal Q = 17,n = 1 model, all in ^K ; the upper triangle of the covariance matrix 
from noise-only Monte Carlo runs, and the lower triangle of the covariance matrix from the best 
fit Monte Carlo runs in /iK . 
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Method 



COBE dataset 
1 year 53x90 
1 year 53+90 
1 year 53 

1 year 53 

2 year 53x90 
2 year 53x90 
2 year NG 

1 year (53A+90A) x (53B+90B) 

2 year 53Ax53B 
2 year 53x90 

2 year (53A+90A)x(53B+90B) 
2 year (53A+90B)x(53B+90A) 
2 year NGAxNGB 
2 year 53+90 



Q? Result 



Reference 

Smoot et al. (1992) 

Wright et al. (1992) 

Smoot et al. (1994) 

Smoot et al. (1994) 

Bennett et al. (1994) 

Bennett et al. (1994) 

this paper 

this paper 

this paper 

this paper 

this paper 

this paper 

this paper 

Gorski et al. (1994) 



Correlation function 
C0BE:cj8 

Genus vs. smoothing 
RMS vs. smoothing 
Correlation function 
Correlation function 
COBE : 1° scale 
Cross power spectrum 
Cross power spectrum 
Cross power spectrum 
Cross power spectrum 
Cross power spectrum 
Cross power spectrum 
Orthonormal functions 
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Y 

Y 
Y 

N 
N 
N 
N 
N 
N 
N 
N 
N 



n 



app 



0.45 
0.65 



n 
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n. 
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n 
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n 
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n 



app 
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1 + 0.23 

■ 1 7+1-3 

■ -L-'_i.i 

. 1 7+0.3 

- 1 c;q+0.49 

■ -'-•'jy-0.55 
1-0.60 
-0.55 



n 
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n„ 



n. 



n. 



app 
app 



n 



app 



n. 
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n. 
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n. 



app 



- i.2i: 

1.15 + 0.2 
: 1.69^ 
: 1.4i: 

: 1.22^ 

: 1.32^ 

: 1.20^ 

: 1.411 

: 1.02 + 0.4 



hO.45 

-0.52 
i-hO.75 

-1.17 
)-H0.42 

-0.46 
,-1-0.39 

-0.45 
1+0.42 

-0.46 
1+0.75 

-1.17 



Table 4: Spectral index determinations 
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Fig. 1. — Power spectrum of the 2 year 53+90 DMR maps (points) compared to the mean it la 
range of Monte Carlo spectra computed for Harrison-Zeldovich skies with an expected Q = 17 /uK. 
The lower band is the iblo" range for noise only Monte Carlos. The lines show the mean power 
spectra for Monte Carlo's with n = —0.5, 0, 0.5, 1, 1.5, 2 & 2.5 all normalized to have the same input 
£ = 4 amplitude as the Q = 17 H-Z case. 
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Fig. 2. — Effective spherical liarmonic index i^jf vs £ for tlie T^ statistics witli a galactic plane 
cut of \b\ > 19.5° (solid curve and points for the binned statistics), and with no galactic plane cut 
(dashed curve). The dotted curve shows i^fj = i- 
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Fig. 3. — Scatter plot showing fitted values for the hexadecupole power and the spectral index n 
using 3 < £ < 19 for 1000 out of 4000 Monte Carlo skies calculated with {Q'rms)^'^ = 1'^ ^^ and 
n = 1. The big circle is fit to the power spectrum of the 2 year 53+90 GHz maps, the big diamond 
is the fit to the mean of the Monte Carlo's, and the dots are the Monte Carlo's. The curves show 
combinations of n and T4 that give Q = 13, 17 &: 22 ^uK (from top to bottom). 
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Fig. 4. — Power spectra normalized to the mean of 17 /iK Harrison-Zeldovich Monte Carlo skies. 
COBE data points from the 2 year NG DMR maps. Models shown as thin curves: n = 1, Q = 17 
/xK is the horizontal line, the best fit n = 1.4 power law is the slanted line, & tilted CDM including 
the effects of gravitational waves with the long dashed curve showing n = 0.96 (predicted by 
0^ chaotic inflation), and the short dashed curve showing n = 0.85 where the tensor and scalar 
quadrupoles are equal (Crittenden et al. 1993). Points with "bent" ends on their error bars are 
from other experiments: FIRS (Ganga et al. 1993), (from left to right) ULISSE (de Bernardis et 
al. 1992), Tenerife (Watson et al. 1992 and Hancock et al. 1994), the South Pole (Schuster et al. 
1993), Saskatoon (Wollack et al.l993), the Python experiment (Dragovan et al. 1994), ARGO (de 
Bernardis et al. 1994), MSAM single subtracted (Cheng et al. 1994), MAX (Gunderson et al. (1993) 
and Meinhold et al. (1993)), MSAM double subtracted. White dish second harmonic (Tucker et 
al. 1993), OVRO (Readhead et al. 1989), OVRO RING (Myers et al. 1993), and the Austraha 
Telescope (Subrahmayan et al. 1993). The open circles above the MSAM points show the effects 
of not removing sources. 



